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We study the low-temperature properties of a one-dimensional spin- 1/2 chain of magnetic impuri- 
ties coupled to a (normal) metal environment by means of anisotropic Kondo exchange. In the case 
of easy-plane anisotropy, we obtain the phase diagram of this system at T = 0. We show that the 
in-plane Kondo coupling destabilizes the Tomonaga-Luttinger phase of the spin-chain, and leads to 
two different phases: i) At strong Kondo coupling, the spins in the chain form Kondo singlets and 
become screened by the metallic environment, and ii) At weak and intermediate Kondo coupling, 
we find a novel dissipative phase characterized by diffusive gapless spin excitations. The two phases 
are separated by a quantum critical point of the Wilson-Fisher universality class with dynamical 
exponent z ~ 2. 

PACS numbers: 75.10.Pq, 75.20.Hr, 75.40.-s 



I. INTRODUCTION 

Magnetic structures of atomic size provide the small- 
est solid state systems in which it is possible to store 
(quantum) information.^ The possibility to build and ma- 
nipulate such atomic-scale magnetic structures has been 
demonstrated in recent experiments using scanning tun- 
neling microscopy (STM),EHS1 a fact that paves the way 
for the realization of spin-devices of nanoscopic size. 

Besides the interest in applications, magnetic sys- 
tems at the nanoscale constitute an excellent play- 
ground to address fundamental questions in condensed 
matter physics. For instance, magnetic impurities in- 
side a metallic host have been found in STM experi- 
ments to display the Kondo effect .'S^ This effect, con- 
sisting in the spin-compensation of a localized magnetic 
moment by conduction electrons in the metal, is one 
of the most paradigmatical phenomena in many-body 
physics.- On the other hand, magnetic atoms inside a 
metal can interact non-locally via the electronic medium 
through the Ruderman-Kittel-Kasuya-Yosida (RKKY) 
exchange coupling,'iSl which is responsible for the mag- 
netic properties of many heavy- fermion materials ""^^ and 
for the giant-magnetoresistance effect in layered mag- 
netic heterostructures.E^ Direct evidence of RKKY in- 
teraction at the atomic scale (i.e., among pairs of mag- 
netic Fe or Co a toms) has been reported recently in 
STM experiments.EE^HlIl Due to the ability to control 
the distance between magnetic atoms using the STM tip, 
the RKKY interaction can be, in principle, tuned from 
ferromagnetic to antiferromagnetic, and this oscillating 
feature has been clearly revealed in recent experiments 
showing spin-polarized STM maps.l^ 

By depositing atoms one by one, STM also enables 
to build magnetic structures where the dimensionality is 
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Figure 1: Scheme of the system. A S* = 1/2 spin chain with is 
coupled to a metallic environment through the Kondo Hamil- 
tonian. The presence of the metallic surface breaks the rota- 
tional and inversion symmetry, inducing in-plane anisotropy. 



gradually changed from the zero-dimensional limit to the 
one-dimensional (ID) case. In particular, linear arrays of 
up to 10 magnetic Mn atoms have been recently built and 
analyzed with STM and inelastic electron tunnel spec- 
troscopy (lETS)"^. Theoretically, low-dimensional mag- 
netic systems are of interest due to the prominent effect 
of quantum fluctuations, which inhibit magnetic order 
and, at low temperatures, lead to quantum phases with 
exotic properties.li^ 

Motivated by the experimental progress described 
above, in this work we study a one-dimensional (ID) 
chain of spin S = 1/2 magnetic atoms (i.e. Kondo im- 
purities) embedded in (or, rather, deposited on) a host 
such like a metallic surface (cf. Fig|T]). The magnetic 
atoms are coupled to each other and to the metallic host 
by means of an anisotropic exchange. The anisotropic 
coupling between the magnetic atoms may be regarded 
as a consequence of, e.g., a Dzyaloshinskii-Moriya inter- 
action resulting from the spin-orbit coupling in the host. 
From a different perspective, this system also provides a 
realization of a ID dissipative system. Indeed, the inter- 
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play between quantum fluctuations and dissipation can 
have important consequences for the phase diagram of 
a system!^ Examples of other physical realizations of 
ID dissipative systems can be found in e.g. resistively 
shunted ID Josephson junctions arrayspSHini Tomonaga- 
Luttinger liquids coupled to dissipative bath j ^SHl^ su per- 
conducting wires coupled to diffusive metal^^^ExUll j^j^^j 
ID ultra-cold atom gases embeded in a Bose-Einstein 
condensate.l^ 

Using a host of analytical methods, which include 
bosonization and renormalization group methods, we in- 
vestigate the effects of the metallic environment on the 
spin chain. At T = 0, we predict that the Tomonaga- 
Luttinger liquid (TLL) phase of the impurity spin chain 
will be destabilized by the (in-plane) Kondo interaction, 
jj^, with the metallic host. Turning on J^, yields a dissi- 
pative phase, whose ground state exhibits long-range or- 
der (LRO) of the in-plane magnetization. Upon increase 
of J^, we predict a quantum phase transition towards 
a disordered Kondo-singlet phase, where the spins of ID 
chain are locally screened by the host electrons, and the 
LRO along the chain is destroyed. 

This work is organized as follows. In Section |ll] we 
present the model for the anisotropic spin chain coupled 
to a normal metal, in Sectionlllllwe show the our main re- 



sults, divided into weak-coupling (Sec. Ill A I and strong- 
coupling (IIIB) treatment of the Kondo Hamiltonian, in 
Section |IV| we present a summary and conclusions, and 
finally. Section IS] is an Appendix. 



II. MODEL 

The Hamiltonian of the spin-chain system embedded 
in a metallic host can be split into three terms: 



from the spin-orbit coupling of the electrons at the sur- 
face of the metallic host. The symmetry conditions for 
the DM interaction to exist are rarely met in the bulk of 
typical metals. However, for an impurity chain that lies 
on a metallic surface where inversion symmetry is bro- 
ken, a DM interaction-term is in principle expected.!^ 
The above spin-chain model, Eq. (|2]), can be mapped 
onto a ID Heisenberg-Ising (XXZ) 5 = 1/2 chain by the 
following transformation.'^ 



e'"^' Sf, 



(3) 
(4) 



where 77 = t'An-^{D/JH) (0 < 77 < 7r/2). Thus, Eq. ^ 
becomes: 



(5) 



/ cos 77 and Jfj — . Although the DM in- 
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teraction in metallic surfaces is smalP^l (i.e., D/Jh <C 1), 
it can dramatically affect the properties of the spin chain. 
In particular, if the initial impurity-chain Hamiltonian 
([2]) is isotropic (i.e., = J^), the above transforma- 
tion (|3| maps it onto an XXZ spin chain with easy-plane 
anisotropy: \Jh/Jh\ = (cos 77) ^ > 1. Therefore, from 
here on, we shall focus on the case of the XXZ spin chain 
with in-plane anisotropy. Under these assumptions, it 
is well known that the low-energy sector of Hamiltonian 
Eq. ^ maps onto the XY— model, whose spectrum is 
described in terms of spi non m odes and exhibits power- 
law magnetic correlations.'^^ISIl (^S^SJ^^^) ~ l^-P"^- 

The coupling between the XXZ chain and the metal is 
described by the following anisotropic Kondo exchange 
interactionP 



H = He +nF + nK, 



(1) 



where He — J2ij a p ^if ^i^j' describes the (static) 
interactions between the magnetic atoms in the chain, 
Hk accounts for the coupling of the chain to the metal- 
lic host, and Tip describes the host electrons. We shall 
not attempt to fully specify the microscopic origin of the 
exchange couplings, J'^^ (ij), since we want to address 
a generic situation in which different contributions (i.e., 
RKKY , direct exchange, dipolar, etc.) may exist. How- 
ever, we shall assume that the dominant exchange cou- 
plings are short-ranged. Thus, we consider the model: 
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Dz ■ (Sj X Sj+i) 



(2) 



= E I^kS^ 4 (R,) ct (R,) - c| (R,) (R,) 



7-L 



S'+4(R,)ct(Rj)+H.c. 



(6) 



where Eqs. ^ and Q have been used. In ([6| every spin 
Sj in the chain is coupled to the fermionic spin-density of 
the bath. The operator cj^ (R) creates an electron with 
spin projection a at position R = (x, y, z) in the metal, 
and for a chain site Rj = (jao,0, 0) with oq the lattice 
parameter of the impurity chain. We also assume an 
anisotropic Kondo interaction with in-plane anisotropy 

Finally, the dynamics of the electrons in the metallic 
host is described in terms of Landau quasi-particles: 



-Hi 



^ek4 (k)c. (k) 

k,(T 



(7) 



For definiteness, we assume here that the index j runs 
along the i— axis (cf. Fig. [ij. The last term in Eq. ^ is 
a Dzyaloshinskii-Moriya (DM) interaction, which results 



where is the electron dispersion and C(j(k) annihilates 
one electron with quantum numbers k and spin projec- 
tion (T. The dots in Eq. ([t]) stand for additional terms. 
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such as spin-obit interactions, whose form needs not be 
specified. 

The physics of model in Eq.([T|) can be understood in 
part by noting that, whereas the Heisenberg interactions 
favor correlations along the spin chain, the Kondo cou- 
pling tends to screen the impurity-spins at every site, 
thus promoting a non-magnetic ground state. However, 
as we will show further below, this intuitive picture where 
the Heisenberg and Kondo interaction compete is not al- 
ways correct. Indeed, the interplay between Heisenberg 
and Kondo interaction is subtle and leads to a counter- 
intuitive cooperation effect in a certain regime of param- 
eters, resulting in a non-trivial phase diagram at T = 
(cf. Fig. |3] below). 

Note that our model significantly differs from previous 
approaches to the ID Heisenberg- Kondo model '^'^ or 
the ID Kondo-lattice model.E3t2Zl Those works assumed 
an entirely ID situation in which the ID spin-chain is 
coupled to a ID electron gas which acts as a host. How- 
ever, our model is closer to the experimental situations 
described in Sec. |T] where the spin chain is embedded in 
or lies on a higher-dimensional metallic host. 

In the following, we focus on the critical properties of 
model Eq. ([T| at low-energies and for different regimes 
of parameters , , and . Throughout we shall 
use units where h — 1. 



transverse direction propagating with velocity u, and 
is parametrized by the dimensionless Luttinger parame- 
^Q^MM\]{ = arccos {-Jh/Jh)] which determines 
the decay of the correlation funtions in the chain, e.g., 
{S+{x)S-{0)) ~ l^r^/^^. The isotropic Heisenberg 
model is recovered for the particular value K = 1/2. 
As discussed above, in this work we focus on the regime 
of easy-plane anisotropy, which corresponds to the con- 
dition that K > 1/2. The ellipsis in Q stands for 
additional perturbations, such as the sine-Gordon term 
^ 008 4$ (x), which are irrelevant in the renormalization- 
group (RG) sense ioi K > ^ and will be neglected. 

The continuum limit of the Kondo Hamiltonian, 
Eq. ([6|, reads: 

-Hk^ f dx its' (x) {x) 



1,3 

dx [e""'"'=S+ (x) s- 



(a;)+H.c.], (11) 



where we have defined qoM = v/'^o and introduced the 
factors of the Fermi momentum fci? in order for the Kondo 
couplings to have dimensions of energy. We have also 
defined the spin-density operator of the host electrons: 



s"(R)^^ct (R) 



1 



Ca' (R) , (12) 



III. RESULTS 

A. Weak coupling scaling analysis 

When the Heisenberg coupling is much larger than 
and J|., it is a good starting point to treat the Kondo 
coupHng, Hk (cf- Eq. [6]), small perturbation to 

the isolated spin-chain Hamiltonian, Hxxz (cf. Eq. [5|. 
In this case, Hamiltonian ([5| can be studied within the 
framework of Abelian bosonization,'^^ which allows to 
map it onto the continuous XY Hamiltonian 



xxz 



1 

2^ 



dx 



- uK [vey 



(8) 



Here 8 [x) , <& {x) are conjugate canonical fields obey- 
ing the usual commutation relations [O {x) , V<i> {x')] = 
inS {x — x'). These fields are continuous in the scale of 
ao, the original lattice spacing in the chain, and are re- 
lated to the original spin operators by 



Sf = aoS"^ (xj) 



2tt 



j-rr/ao 



cos 2$ (xj) 



(9) 



ao. 



S] = aoS' (x,) = -^V$(.T,) 



COS 2$ (xj) , 
(10) 



where Xj = juo is the position of the j-th spin. 
The model (Isl describes ID spinon excitations in the 



with cr" (a = x, y, z) the Pauli matrices. From Eqs. ( 10 1 



and (111 we note that the operator V$ [x) couples to the 
uniform component of the spin-density in the electron- 
gas [x) = s' (x), and the operator cos 2$ (x) couples 
to the staggered component sf (a;) = e^^^/°-ogZ 
the other hand, the operator e"*®^^^-' couples to the stag- 
gered magnetization sj (x) = g^^^/i^og- (^j;). These con- 



tributions to Eq. (Ill have different scaling dimensions, 



and we treat them independently in the following analy- 
sis. 

Next, we assess the stability of the Tomonaga- 
Luttinger (TLL)( phase, which is described by Hamilto- 
nian in Eq. (|8| . To this end, we have considered the lead- 
ing order corrections to the free energy per unit length of 
impurity spin-chain. The technical details of this calcu- 
lation can be found in Appendix |^ We shall not pursue 
the stability analysis beyond the leading order in this 
work, as our focus here is on the phase diagram in the 
K>l/2 (i.e. Jjj > Jfj) regime, for which, as the follow- 
ing discussion demonstrates, there is only one relevant 
Kondo coupling, namely jj^. A more complete analy- 
sis will be reported elsewhere To leading order in the 
Kondo couplings, for temperatures T <^ J^, we find: 



9l,u 



+ 



,2 

:} Z.S 



267r2u2fc^ 



24fcj.a27r3 \ u J 



rp6 



(^Y'^ J.1+2K 



+ 9l,s 



/7rao\i/2^^i+i/2X 



25fcFa27r3 
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(13) 
(14) 

(15) 
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where L — +00 is the impurity chain length and the di- 
mensionless couplings g^^u = Jx/vEkF, gz,s = JK/vphp 
and g^^s = Jii/'VFkFj have been introduced; , and 
are non-universal numerical coefficients. The stabil- 
ity of the TLL phase with respect to the perturbation 
T-Lk can be now assessed by comparing the scaling with 
temperature of A^F and the free energy of the spin-chain 
chain, fp, when decoupled from the environment, which 
behaves aPI^Fo - T^. Thus we look for divergences in 
the perturbative corrections to AF/Fq as the tempera- 
ture is gradually decreased towards the ground state (i.e. 



T = 0). From (151, it can be seen that the term oc gj^ ^ 



yields a divergent contribution to AF/Fq, which signals 
an instability of the TLL phase. 

To make contact with the renormalization group 
(RG), we shall define the effective couplings gzu{^) = 
7.,ne-(i/2)^, ^^^^(^) ^ 5.,.e(i/2-i^)^ and (£) ee 



g±,se 



{l/2-l/4:K)£ 



, where £ = \n ( Aq /T) and Aq 



is the 



high-energy cutoff for the effective low-energy description 
of the spin-chain in terms of Eq. ([s]) . Decreasing the tem- 
perature a bit towards the ground state (i.e. T = 0) can 
regarded as an infinitesimal change of £ ^ £ + S£, and 
the corresponding change ('fiow') of the effective cou- 
plings can described by the following set of differential 
equations: 



dgz,u (^) _ 1 

d£ ~ 2^"'" ^ ' ' 



dg±,s {£) 
d£ 



1 



1 



1 
2K 



g±,s (e) ■ 



(16) 
(17) 
(18) 



Alternatively, we can regard these equations as describ- 
ing the change in effective (dimensionless) couplings of an 
equivalent (coarse-grained) system with a reduced high- 
energy cutoff A{£) — e^^Ag. This interpretation means 
that for in-plane anisotropy where K > 1/2, the the 
couplings gz^u (^) and gz^s (i) decrease as the system is 
coarse-grained by integrating out the high-energy degrees 
of freedom and become irrelevant (in the RG sense). In 
other words, the terms in proportional to those cou- 
plings yield subleading corrections and therefore can be 
neglected as T — >■ 0. On the other hand, g±^s (^) is a rel- 
evant (in the RG sense) perturbation, which, as T 
yields an dominant correction and destabilizes the TLL 
phase of the spin-chain described by Eq. ([s]) . Note that 
both gz.s {£) and g±^s {£) have the same scaling dimen- 
sion at the Heisenberg point {K = 1/2), where they 
are marginally relevant, and a higher order perturbative 
analysis is required to fully assess the stabilty of the TLL 
phase. Bfil 

Finally, as explained in the Appendix [X] in the deriva- 
tion of Eq. ( 15 1 we have assumed that the host is a nor- 
mal metal. This means that it exhibits an ohmic spec- 
trum of magnetic excitations over a broad range of mo- 
mentum transfer along the spin-chain direction. Alter- 
natively, this situation can be mimicked by an infinite 



set of fermionic baths, each bath being locally coupled to 
only one impurity spin (cf . Fig. |2| , which yields a local 
(i.e. momentum independent) ohmic spin response: e.g. 

{x,t),x^^{x,t) ~ 6{x)/t'^. As it will be discussed in 
the next section, this model allows us to treat the rele- 
vant Kondo coupling jj^ in a non-perturbative way. In 
particular, it captures the important (non-perturbative) 
feature that the magnetic moment of the impurities will 
be fully screened by the metallic environment at large 
Jj^ ^ J^. In the above analysis, the need for a non- 
perturbative treatment is evidenced by the fact that even 
an infinitesimal value of will destabilize the TLL 
phase for in-plane anisotropy {K > 1/2). However, dif- 
ferent from the single-impurity Kondo problem, we will 
see below that the RG fiow does not proceed from the 
TLL phase into a strong coupling Kondo-screened phase 
in a straightforward manner, but rather, another phase 
of dissipative nature intervenes between the TLL and the 
Kondo phase. 



B. Strong coupling analysis 

1, Derivation of an effective ID model 

As mentioned before, when jj^ fiows to strong cou- 
pling, the perturbative RG approach used in the previous 
Section III A is no longer valid, and we need to study the 
physical properties of the spin chain in a different way. 
The approach used in this section is motivated by the fol- 
lowing arguments: ii) the analysis made in the previous 
section indicates that the most relevant coupling of the 
spin-chain to the host stems from a local kind of coupling 
to the metallic host and, ii) even at strong-coupling, for 
a two- or three-dimensional host, the interference of two 
Kondo screening clouds belonging to spins located at a 
distance d decays rapidily when d is of the order of a few 
Fermi wavelengths'^^"'*'* (i.e. ~ fc^^)- Experimentally, 
this is confirmed by the behavior of the STS Fano line 
shapes in magnetic Co atoms deposited on Cu(lOO) and 
separated by distances d > lOA, which are identical to 
the single-impurity STS line shapes .^^^ 

These arguments suggest to perform an approxima- 
tion on Hamiltonian T-Lf, assuming a set of independent 
fermionic baths (i.e., semi-infinite ID chains, cf. Fig. [2]) 
coupled to each single spin in the chain 



H.C. 



(19) 



Here, the index runs along the i (y)-axis. Eq. (19 1 
is a valid description in the limit kpao 3> 1, consistent 
with the local limit of Sec. IIII Al Note that in this local- 
bath approximation, since the fermionic chains are not 
connected between them, the indirect RKKY exchange 
interaction mediated by conduction electrons is not dy- 
namically generated at second order in the Kondo cou- 
pling and it is necessary to incorporate it explicitly in the 
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Figure 2: Schematic diagram of the spin-chain coupled to 
the metallic host, after the independent-bath approximation. 
Due to the locality of the dynamic spin-susceptibility, which 
decays on distances of the order of ~ fc]^^, the interference of 
Kondo-screening clouds is negligible in the limit kpao 3> 1. 



Heisenberg exchange parameters jj^,Jfj. In the follow- 
ing, we assume that both J^,J^ already contain these 
corrections. 

The advantage of the independent-bath approximation 
Eq. ( 19 1 is that it allows to use powerful analytical meth- 



ods which have been applied successfully to describe the 
single Kondo-impurity problem. In the following, we im- 
plement the Abe lian bosonization approach to the Kondo 
problem.liSliiniHiZl Xo avoid confusion with the previous 
Sec. |III A| note that here bosonization is implemented to 
describe the fermionic ID chains, and not the spin chain. 
At low energies the Hamiltonians H f and Hk become in 
the bosonic representation^^ 



/oo 

(0) 



(20) 



n 



K 



g-J\/20f^(O) pi\/205(O) 



'0 



27r6o 



C {y) , (y) 



(21) 



obey 



where the chiral fields 
the commutation relations 

insignijj — y') Si,jS^^rii and are related to charge 
and spin density-fluctuations through the rela- 

tions p^y) = E.RlAy)R^,Ay) - -^v</)ff,(y) 

and m, {y) = ^4. iv) (?/) = -^^^t 
respectively.^^ In Eq. (20 1 vp is the Fermi velocity, and 
in Eq. (21 1 6s = ta.ir'^'^{TTpQj^bQ/A) is the scattering 



phase-shift associated with the potential 
po = {2TrvF)~^ the conduction electron density of states 
at the Fermi energy, and 69 the lattice parameter in the 
fermionic chain. For simplicity we assume these param- 
eters to be identical for all chains. We then introduce 
the (Emery-Kivelson) unitary transformation^ 



U — exp 



-J7 



E^»''^5(o) 



(22) 



under which the bosonic field V0f^ (y) and the spin op- 



erator transform as 

Z^t (y) U = [V05 {y) + 5 (y) 27r^S:] , (23) 
U^S+U = S+e'^*^^^(°K (24) 

Upon this transformation, the model Hamiltonian, 
Eq. (g transforms as ?i = WnU = Hc + Hk + Uxxz, 
with 



T~Lc — T~Lc, 

^ TTpo V2tt 



-»(y2-7X,(0) 

St — +H.C. 



27r6o 



(25) 



, (26) 



iixxz ~ E JhSiSl^i' 



J2 ■Ife'A'Ptio)-^?+..sW]s+S-,,+U.c. 



(27) 



where we have defined 5g = Sg — nj/2\/2. Note 
that in the transformed representation, the quantum 
dynamics of the bath [represented by the chiral field 
(j)^g (0)1 appears explicitly in the Heisenberg term ^ 
(^S+Sr^^e'^l'l't(o)~'l'f+i.Ao)] +h.c.). Physically, this 

means that the Heisenberg interaction is now "dressed" by 
the spin-density fluctuations of the electron gas. Note 
that the independent-bath model ( 19 1 is crucial to im- 



plement bosonization along the chains, and to put these 
ideas on a clear mathematical framework. 

Up to now the parameter 7 appearing in Eq. ( 22 1 
remains completely arbitrary. We now set 7 = \/2 in 
Eqs. (26 1 and (27 1, and the tranformed Hamiltonians 



simplify to 



2Ss .,V0^ 
Trpo V2n 



(0) 



7-L 

27T ' 



(28) 



^ :ke»V2[C(o)-0f+i,.(o)]5+5-^^ + H.c, (29) 

where now T-Lk is equivalent to the spin-boson model with 
Ohmic dissipation with 5s related to the dissipative 
parameter a in the context of macroscopic quantum co- 
herence through a = (^Ss/t^ , and with the in-plane 
Kondo interaction playing the role of a magnetic field 
along the x— axis ~ — J^/27r. That model describes 
a quantum phase transition from a phase with a "frozen" 
spin state (either S- = +1/2 or = -1/2) for a > 1, to 
a phase with an "untrapped" spin state for a < 1, where 
spin-fiips induced by Sf proliferate. The strong-coupling 
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regime of the single-impurity Kondo Hamiltonian there- 
fore corresponds to this last case, where (not too close to 
the transition) the Kondo temperature i^Mmmi 



VF 



a/{l-a) 



(30) 



The special case a = (i.e., 5s = 0) was analyzed in 
Ref. [47j and represents a particular limit where Hk can 
be diagonalized in the eigenbasis |±)j of the operator S'f . 
Note that the condition = implies that the origi- 
nal phase-shift is 5s = 7r/2, corresponding to the unitary 
limit J^b/vp — )• oo. Unfortunately, the unitary limit is 
not consistent with the well-known local Fermi-liquid de- 
scription of the strong-coupling Kondo fixed point .'^In- 
tuitively, in the limit 5s — the coupling to the bath van- 
ishes and Hk reduces to a Zeeman Hamiltonian, which 
is not equivalent to the Kondo problem. However, as 
shown by Kotliar and Si, a physically correct description 
of the strong-coupling limit is recovered by performing 
second-order perturbation expansion in ^s.!^ Physically, 
this is equivalent to reintroducing the coupling to the 
bath. 

In the case of plane-anisotropy, the Kondo couplings 
satisfy > J^. This implies that in order to perform 
an expansion around the point Sg — 0, the microscopic 
parameters of the model should be in a regime such that 
strictly speaking the use of bosonization is not justified 
(i.e., the interactions are of the order or bigger than the 
Fermi energy). However, since this approach has been 
shown to successfully capture qualitatively the main fea- 
tures of the strong-coupling Kondo fixed point we ex- 
pect our approach to be correct only at a qualitative level. 
In addition, since 5s ~ implies a ~ 0, we will use Eq. 



(301 to identify Tk 



in what follows. 



An effective low-energy Hamiltonian valid near 5s = 0, 
and for the case where the Kondo interaction dominates 
over Heisenberg exchange, i.e., ^ {'^H''^h} '^^^ t)e 
obtained expanding Eq. ( 28 1 at order 5s and projecting 



onto the lowest energy level on each site 



n 



I 

K = 



1 27r 



(V</.5(0))', (31) 



where we have introduced the projector operator on the 
subspace spanned by |-),^ , i.e., P_ = Hi (|-), and 
where we have neglected a constant term J^/27r. The 
effective magnetic field opens a gap of size A — 2hx 
in the spin excitation spectrum and consequently the 
spin degrees of freedom are "frozen" in the lowest energy 
configuration |— )j. In contrast, spin-density fluctuations 
in the bath remain gapless and their dynamics becomes 
dominant at low energies. Projecting Hxxz onto this 
basis yields 



xxz 



= P-U 



xxz 



p_ 



2 42 



cos 



V2[</.5(O)-0f+i^, (0)] 



(32) 



Heisenberg term 
effective interac- 
encoded in the 
_ ]^ 3 (0)] . This is an im- 
work, complementary to the 



In this representation, the 
(Sj^S'j^j -I- h.c.) induces an 
tion between neighboring baths 
term- cos ^/2 [(/.^ (0) - 
portant result in our 
situation analyzed in Sec. |III A[ where the opposite 
limit ^ was studied. In that case, the bath 
was integrated out, and we studied the stability of the 
Tomonaga-Luttinger liquid phase to leading oder in 
perturbation theory. Here, we do just the opposite: we 
keep the degrees of freedom of the bath and eliminate 
the spin degrees of freedom. 

Although a general derivation of an effective low- 
energy model (i.e., for an arbitrary value of 7), is beyond 
the scope of the present work, it is worth noting that the 
effective coupling between fermionic baths that appears 



in ( 32 1 is a physical feature that goes beyond our particu- 
lar derivation for 7 = v^. This can be understood using, 
for example, similar arguments as those leading to the 
Nozii'^ceres' local Fermi-liquid. Indeed, when ^ J^, 
a natural approach is to start from the Kondo singlets 
at neighboring sites i and i -I- 1, i.e., |G,;) and |Gi+i) re- 



h.c. 



spectively, where|Gi) = c\^^^j - c\^.^ 

The perturbation %' = -< I]i={i,i+i} (ca,crC/o,(T 
acting on these neighboring singlets produces virtual 
excitations to the n;o = 1 triplet subspace, at order 
{t/J^) on each one of them. Eventually, the Heisen- 
berg interaction {Si^ +h.c.) restores the initial 
singlet ground states \Gi) and |Gi+i) and, as a net re- 
sult, virtual processes generate an effective spin inter- 
action — (t/ J^^^ (^cl ^ ^Ci^i^i^i + }i.c^ between the 

second sites in the chains i and i + 1. Bosonizing this 
induced effective interaction yields a term of th e fo rm 
- cos V2 (0) - , (0)] , analogous to Eq. 

Next we can proceed to derive an effective ID model, 
which first proceeds by the integrating out of the modes 
(/)f^^ (y) for y ^ 0. This can be done exactly using the 
functional integral representation of the partition func- 
tion, and generates a term of the form ~ \u}m\ hi the ef- 
fective action, which stems from the (ohmic) dissipation 
induced by the coupling to local bath (cf. Refs. I16|41[ 
for details). The resulting Euclidean action of the sys- 
tem reads: 



S' — S'q + S'ff, 



(33) 



^,JdrV^ {r)Y 



2TTEn 



rP j± 

S'h^YIj^ dr ^ cos [(f^ (t) - ip, (r)] 



(34) 
(35) 



where we have defined the more compact notation for lo- 
cal field ipi = ^/2(f>f ,, (y = 0), and where we have used the 
equation of motion of chiral fields, i.e., dr^i — ivp^I Lpi = 
to express V(y9i in terms of drfi- In addition, in (34 1 we 
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have defined the parameter 



Ec 



2~5V 



suit (cf. Eq. B14l C+~ (0,t) r~2, as expected fo r the 
local Fermi-liquid description of the Kondo problemP^I^ 
This slow decay is a consequence of the ohmic dissipation 
(36) term ^ |w„i| in Eq. (34), inherited from the dynamics of 
the semi-infinite fermionic chain. This behavior is con- 
sistent with the strong-coupling Kondo picture where, at 
long times, the spin degrees of freedom are merged with 
those of the fermion bath.**^ On the other hand, the expo- 



where the singularity at i5s = is a consequence of the 
unphysical unitary limit mentioned above. 

Note that the effective model Eq. (33 1 is formally nentially decaying non-local correlation [cf. Eq. (B15l] 



equivalent to the action of a ID Josephson-j unction ar- 
ray with local Ohmic dissipation, with the phase 
of the superconducting order parameter at site i, jj^ 
the Josephson coupling^ and Eq the charging energy 
with respect to the ground .'^Sl It can be also brought 
to a form equivalent to a ID 0(2) dissipative quan- 
tum rotor model if we w rite it in terms of N,; (r) = 
(cos Lpi (r) , sin tp^ (•^)) | 25 | 55| l57] rpj^g ^gj^^ ^Yig^^ Hamiltonian 
Eq. ([T|) can be mapped (in the limit J^, ^ J^, J^) 
to these dissipative models is an important result of our 
work which shows interesting underlying connections be- 
tween appartently different physical situations. 

To appreciate the physical consequences of the effective 
model in Eq. ( 33 1 , we concentrate on the transverse spin 
correlation function C+- In.r) = (r^5^;„ (r) S"," (0)). 
In the transformed representation, this correlation 
is evaluated as (cf. Appendix IbI) [n, r) — 



Y Near the strong-coupling single- 



impurity Kondo limit 3> { , Jfj } , we obtain the re- 



1 g-n/4c 

C+" {n,0) = -—-., 

2 n + 1 



(37) 



with = 1/ln the correlation length, indicates 

that the spins are not spatially correlated beyond a dis- 
tance supporting the idea that in this limit the spin 
chain realizes a disordered phase of nearly independent 
Kondo singlets. 



Hubbard- Stratonovich decoupling and RG analysis of the 
effective ip^ — theory 



The properties and phases of action ( 33 1 can be inves 



tigated introducing an auxiliary bosonic field ipi (r) to 
decouple the Heisenberg term (cf. e.g. Refs. 
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E 



dr cos [ipi (r) - (pi+i (r)] 



dr J2 (r) [J-'] (r) - ^ [v^* (r) e^'^^(^) + A (r) e-''^'^^) 



(38) 



where we have defined the matrix [Jl-- — 
^Jji Then, the partition function 

reads Z = Zq^V[iP] e-'^'^l , where 

I 1 r^/^ r ' 

- In / exp - ^ / dT (r) e'^^^^^ + h.c. 

(3°9) 

is the effective action for the auxiliary field tpi {t). Here, 
the notation (. . . )q means average with respect to the 
local action ( 34 1 . A cumulant expansion of the last term 



in Eq. (|39|) to order xp'^ yields 



2N, 



4! 



E 



l/T 



dr lip, (t)| , 
(40) 



where the compact notation q = (fc,Wm), with uj„i = 
2i:mT the bosonic Matsubara frequencies,'^ has been 
used, and where Ni is the number of spins. Here we 
have defined the Gaussian propagator 



16 



-27 



H 



El 



where 



32 7re~T^ 
c 



(41) 



I dr e'"-- /T,e'*"("i)e-''*"("^'\ (42) 

(43) 
(44) 



with 7b — 0.577 . . . the Euler gamma constant and c = 
21.8 ... At the mean-field level (i.e., the saddle-point ap- 
proximation to Eq. 40 I , this model describes a QPT when 
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Disordered Kondo Phase 



Dissipative Phase 



K^=i/2: 



Luttin^ei^l^id 



- Jt 



Figure 3: Schematic T — phase diagram of the dissipa- 
tive ID Kondo-Heisenberg spin chain. The blue line corre- 
sponds to the Tomonaga-Luttinger liquid phase, which is un- 
stable against a small perturbation J^. The dashed line is the 
critical line separating a phase of disordered Kondo-singlets 
(upper white area), from a dissipative phase (shaded bottom 
area), characterized by gapless diffusive spin excitations. In 
this phase, at T = 0, the dissipative dynamics induced by 
the metallic environment stabilizes long-range order in the 
transverse magnetization. 



reaches the critical value jjj ^ = QAEq/ttc^^'^ cx jj^. 
In Fig. |3] we show schematically this critical line as a 
dashed black line, separating the disordered phase with 
(i/'i) = 0, corresponding to the disordered Kondo phase 
described in Section IIIB [cf. Eq. (37l], from the or- 
dered phase with (?/;,;) ^ 0, corresponding to the dissipa- 



tive phase of Section III A Physically, when > ^ 



the Heisenberg interaction is large enough to induce long- 
range coherence in the transverse magnetization (5*^^) oc 
(^e'"^') oc {ipi) along the spin-chain. On the other hand, 
when < ^ the "charging" term induces large 
quantum fluctuations of the field Lpi (i.e., local spin-flips 
induced by the term J^) and therefore tends to destroy 
the ordered state. In this representation, the competition 
between the couplings jj^ and is transparent. 

Beyond the mean-field level, the quantum 
critical properties of this model, generalized 
to describe a TV— component field {-^q (a;)} — 
{ipi (x) , "02 {x) , ■ . ■ , '4'N {x)} in d dimensions, are 
well-known and have been studied in the context of an- 
tiferromagnetic instabilities of Fermi-liquids,^^ ^■^ using 
the framework of the Hertz-Moriya-Millis theory.ESlMI 
This theory describe quantum fiuctuations of the order 
parameter {x), with critical dynamical exponent 
z — 2, around the Gaussian fixed-point [i.e., m = in 
Eq. (40 1]. A standard momentum-shell RG procedure, 
performing a two-loop expansion in u, and in the small 
parameter e — A — D, where D — d + z, allows to obtain 
the RG-fiow equations for the parameters of model 



(40l,L 



dr (£) 

dt 
du{t) 

de 



2r {£) + u {t) 
eu {€) - {£) 



N + 2 A 



D 



6 

N - 



1 



-r(£)' 
Ad 



6 {l + r{i)f 



(45) 
(46) 



where Ad = 27r(^/2)/r (d/2) is the surface area of the 
£)— dimensional sphere, and where the units are such 
that the high-energy cutoff of the theory ( [40| is A = 1. 
For e small and positive {d < 2), this RG-fiow is con- 
trolled by a Wilson-Fisher fixed-point located at r* — 
-e {N + 2) /2 {N + 8) ,u* = 6e/AD {N + 8). 

The stability of the Wilson-Fisher fixed point can be 
studied upon expansion of Eqs. (45 1 and (|46| in the small 
deviations dr = r — r* , Su = u — u* , and allows to obtain 
the eigen-coupling equation dwi/di = XiWi, {i = 1,2), 
with eigenvalues Ai = 2 - e (A^ + 2) / (A^ + 8) + O (e^) 
and X2 = —e + O (e^) , determing the critical exponents 
of the transition .^S21£2) Although for d = 1 the parameter 
e = 4 — (d + z) = 1, which is not small, recent quan- 
tum Monte Carlo simulation^^Zl have shown evidence of 
a QPT in Ai^-spin chains subject to local Ohmic dissi- 
pation, with critical exponents in good agreem ent with 
those predicted using the the e-expansion .123221 

In the ordered phase (ipi) ^ ocurring for > J^^,, 
since the U(l) symmetry of Eq. (39 I is spontaneously 
broken, a Goldstone-mode arises from smooth fiuctua- 
tions of the phase of the order parameter -0^ ~ e"*', 
and it becomes necessary to check the stability of the 
ordered-phase. To that end, we return to Eq. (40 1 and 



perform an expansion in small fiuctuations of the phase 
bdi = 'di — 'do, around an arbitrary value ??o- At Gaussian 
order in S'di we obtain the effective action 



G-J (q) 



1 1 

2M. 



E 

q 



G-J{ci)\d{cDY 



El 



1602 



(47) 
(48) 



where G~g (q) is the propagator of the Goldstone mode, 
and where ipo is the saddle-point solution of ( [40| . This 
propagator describes a gapless phase characterized by 
diffusive (z = 2) excitations of the field i?(q), and by 
correlation functions {x,t) = {TrS~^ {x,t) S~ (0)) 
decaying as C"' (2^,0) ^ ^ at long dista nces , and 
C+- (r) - |rp^/^ at long times (cf. Appendix B2l. Us- 
ing Eq. ( 47 1 to evaluate the average of the order param- 

eter (V;,) = (e^^O = V-oe""^"'' 



we obtain the result 



(V'i) = "00 exp 



(for T = 0). 

(49) 



Interestingly, due to the presence of the dissipative term 
^ |wm|, the (Gaussian) fiuctuations of the spin-chain 
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are strongly suppressed relative to the isolated (XY) 
chain. Indeed, contrary to the case of isolated ID sys- 
tems, quantum fluctuations do not destroy the LRO be- 
cause the effective dimensionality of the quantum system 
is D — d + z — 3, larger than the critical dimension 
Dc — 2 determined by the Gaussian theory for the Gold- 
stone mode Eq. (47 1. Note that this is not in contradic- 
tion with the Mermin- Wagner theorem,''^ which predicts 
the destruction of LRO at T = in ID systems with short 
range interactions. In our case, due to the presence of a 
higher dimensional fermionic bath which induces long- 
ranged (imaginary) time correlations, the system cannot 
be considered strictly one dimensional. Therefore, the 
Mermin- Wagner theorem is not applicable to our system. 

At this point, it is interesting to note the connection 
with the weak-coupling TLL description of Sec. |III A| 
which becomes apparent using, for instance, the self- 
consistent harmonic approximation (SCHA) method. 
This method consists in finding the optimal propaga- 
tor Gj^j^g^j (q) of a trial Gaussian action, such that the 
variational free-energy of the system is minimized. For 
the case of Sec. 



turbation to the isolated spin-chain. In that case, we 



III A 



when K > 1/2 the only rele- 
vant variable is the field 6 {x), describing the transverse 
spin-excitations of the spin-chain weakly coupled to the 
metal, and consequently the trial action writes compactly 
as S'triai = 2irEq<^t~riai(q)|0(q)r Here we do not 
show the derivation of the SCHA equations, and refer the 
reader to Refs. |25I27] , where this method was applied to 
closely related ID TLL systems in contact to dissipative 
baths. The optimal propagator found in those works is 

Atrial (q) = V \^m\ + K [up + U?Ju) /tT, with Ty - J^. 

The similar forms (in the limit q — >■ 0) of Gj^j'^^j (q) and 
(q) in Eq. ( 48 1 suggests that the dissipative gap- 
less phase with z = 2, obtained in the strong-coupling 
regime, is also stable in the weak-coupling regime (see 
lower part of Fig. [s]). The thick solid line at the bot- 
tom (i.e., = 0) corresponds to the gapless TLL phase 
described by the Hamiltonian ([8|. As shown by the RG- 
fiow Eq. ( |18| , this line is unstable against a vanishingly 
small perturbation J^. 



IV. SUMMARY AND CONCLUSIONS 

We have studied the quantum critical properties of a 
S" = 1/2 spin chain described by the anisotropic XXZ 
Hamiltonian ([s]), coupled to a metallic environment via 
the anisotropic Kondo model ([6]). We study this system 
in different regimes of parameters using various analytical 
approaches (i.e., Abelian bosonization, renormalization- 
group method, analysis of the Ginzburg-Landau func- 
tional, etc.), and obtain the quantum phases at T = 0. 
We focus on the case of plane-anisotropy in Hxxz and 
Hk, that favors the couplings in the transverse direction. 
In real systems, this kind of anisotropy is expected in a 
system with lack of inversion symmetry (see Fig. [T|) . 

First, we investigated the weak-coupling regime <C 
(see Sec. Ill A I, where Hk can be considered a per- 



derive a set of RG-flow equations [cf. Eqs. (16)-(18l] at 
first order in in the parameters and . This RG-flow 
is determined by the value of the Luttinger parameter 
K of the chain, and is drastically different the RG-flow 
expected for the single-impurity Kondo problerrP. Far 
away from the isotropic SU(2) point K = 1/2, while 
becomes an irrelevant coupling (in the RG sense) , be- 
comes relevant and destabilizes the Tomonaga-Luttinger 
liqued flxed-point (cf. thick bottom line in Fig. [3|. This 
situation is analogous to other Lutti nger- liquid systems 
coupled to a dissipative environment .12S123 

In the opposite limit ^ J^, we consider the prob- 
lem starting from the limit of decoupled Kondo-screened 
impurities. We show that the effective model for coupled 
Kondo impurities [cf. Eq. (33l] is formally equivalent 
to that of a ID Josephson-junction array with on-site 
Ohmic dissipation, which is known to undergo a quan- 
tum phase transition as a function of the dissipation 
parameter.l^SI To extract the properties of the spin chain 
in this limit, we derive an effective Ginzburg-Landau the- 
ory [cf. Eq. (40 1[ in terms of a U(l) bosonic Hubbard- 



Stratonovich fleld -0^, and stud y the critical properties 
around the Gaussian flxed-point .1^21621 Physically, the or- 
der parameter ipi describes coherence in the transverse 
magnetization {S^') cx {ipi) along the direction of the 
spin chain. In the limit 3> J^, the system is charac- 
terized by a vanishing order parameter (-0^) = 0, which 
we interpret as a manifestation that the system realizes 
a disordered Kondo-singlet phase, where the spins are 
screened by their local fermionic bath. Upon the in- 
crease of , the system experiences a QPT (dashed line 
in Fig. [3]) towards a phase with LRO, characterized by a 
non-vanishing order-parameter (tpi) ^ at T — 0. Inter- 
estingly, the dynamics of the emerging Goldstone mode 



[cf. Eq. (47l[ is not able to destroy the mean-fleld solu- 



tion, which would be the usual situation for ID systems. 
This anomaly is because of the dissipative character of 
the dynamics of this mode. 

Our results, summarized in the T — phase-diagram 
in Fig. [3] crucially depend on the assumption of in-plane 
anisotropy in the problem. As shown in Refs. [30 31], 
this assumption is physically reasonable due to the pres- 
ence of Dzyaloshinskii-Moriya interactions (induced by 
spin-orbit coupling at the surface of metals) which break 
the SU(2) invariance, and enhance the transverse Heisen- 
berg coupling [cf. Eq. Other physical mecha- 

nisms producing other types of anisotropy, such as mag- 
netocrystalline effects in the metallic host are beyond the 
scope of the present work. 

An important approximation we used is the assump- 
tion of locality in the spin-response of the metal. The 
conditions when this approximation is justifled are dis- 
cussed in detail in Sec lIII Al 

It would be very interesting to verify our predictions on 
experimental level. Experimentally, the phase diagram 
could be studied with STM techniques either by varying 
the strength of the Kondo exchange coupling (i.e., by 
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growing the ID spin chain on the top of a decoupling 
layerl^, or by changing the distance between spins, which 
has the effect of changing the magnitude and sign of the 
exchange interaction. 
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where Zq = Tr e"^"/^ is the partition function of the 
uncoupled host-chain system, and Ho = "Hxxz +Hf (cf- 
Eqs. [8]and[7|, and t is the Matsubara imaginary time.l^ 
The average {A)^^ = Tr e~'^°/'^A/TT e"^"/'^ stands for 
the thermodynamic average of an operator A over the 
Gibbsian sensemble defined by Hq. 



Appendix A: Leading corrections to the spin-chain 
free energy 

We begin by considering the partition function 
Z = Tr e-(«o+«K-)/T ^ (r^e" -^o'" drn^<-)) , (Al) 



The free-energy of the system (removing the bath 
contribution) is given by F — i^o = —T\n{Z/Zo) = 
-Tln^T^e" . Upon expanding the expo- 

nential to the lowest order non-trivial, we obtain the lead- 
ing free energy correction: 



AF ^ -Tin 



2! 



l/T 



dTidT2 {TrUK (n) Uk (t2))o + 



dYidY2 



{TrS^ {v,)S^- {v2))oX'' (ri2) 
(A2) 



where we have introduced the compact notation r = 
{x,t), and ri2 = ri — r2. We have also defined the 
host spin response function, x"'' (R-i ~ R-2i '''i ~ '''2) = 
{T^s- (Ri,Ti)s'' (R2,T2))q - X^io(Ri - R2)<5(Ti - ra), 
X°L(R-i — R-2) being the static part of the spin-response 
of the host, which has been already included in the static 
spin-chain exchange couplings Jji and Jfj that deter- 
mined K . In deriving the above expression, we have used 
that the U(l)xZ2 symmetry of the XXZ chain implies 
that {rr[S+{vi)S+{v2)])o = {Tr[S-{v^)S-{v2)])o = 0, 
etc. Thus, even if the spin-response of the metallic host 



lacks of the in-plane spin-rotation [U(l)] symmetry of 
the spin-chain due to e.g. spin-orbit interactions (i.e. 
X^^(R, t) ^ 0, et cetera), the leading corrections to the 
free energy are insensitive to it. However, one factor that 
complicates the weak-coupling analysis is the fact that 
the spin correlation functions of the host electrons have 
very different behavior as compared to the spin correla- 
tors of the uncoupled spin chain. In particular, the con- 
formal invariance (e.g., homogeneous scaling upon space 
and time rescahng) present in the XY model, Eq. (Is])), 
implies the spin correlators of the uncoupled chain readP^l 
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. u J 



h.c. 



(w) 



(r.5+ (r)5- (0)) = 



( 7TX \ 
aa J 



«^ sinh^ ( -^(-+— ) ) sinh^ 



2^0 sinh* 



L ^ ^T{x+^ur) ^^ ^^^^^ (^^^^ 



IV 



(A3) 
(A4) 



where we have kept the leading terms at large dis- 
tances and times. The uniform component of 5^ 
^-ie{x) j,Qg 2$ (a;)] is less relevant than the staggered part 
gia:7r/aog-ie(a:j)j ^^^-^j ^j-^-^ ncglcctcd in what foUows. 
Estimating the temperature of the various contribu- 



tions to Eq. (A2) is possible by realizing that the func- 



tions that multiply the staggered (oc cos(7ra;/ao)) part of 
the spin-chain correlation functions vary very in space 
slowly compared to the oscillatory factors e*^'"*'^*'^/'^'''^. 
For instance, this allow us to approximate the integral: 



/ 



d:rid:r2e'('°-±^/'^°)(^i-^^)F^(a:i - Xz, ri)x+- (xi - X2,n - r^) ~ LFk{Q,t^ - T2)x^-{qDM ± -,ti - T2) (A5) 



where 
Fk{x,t) 



sinh 



(A6) 



X 



[Qtt) = j dx 



(x,t) 



(A7) 



The spin susceptibility of the host, x (Qi ''') is a rather 
complex object. However, for a normal Fermi liquid, the 
long time dynamics is dominated by particle-hole exci- 
tations, which yield an ohmic behavior (cx u, where lo 
is the small excitation frequency) for a wide range of 
momentum transfer Q < 2kF- Thus, after Fourier trans- 
formation. 



x+-(0,T»c.-i)^G+-(g) 



T2 



(sin -kTtY 



(A8) 



where lOc ^ Ep, that is, the Fermi energy of the host. 
The contribution of the uniform part of the — Kondo 
coupling can be treated in a similar fashion, with the 
difference that the spin-chain fluctuations in this case 
couple to the spin response of the host at Q = 0. 
Using these approximations in Eq. (|A2[)) yields 
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rp6 



+ 5i 



/7moy/2^yi+i/2K 



\ u J 



(A9) 



where we h.£ive defined the diinensioiiless couplings Qz^u — 
J^/vphp, gz,s = Jx/vEkp and g^^g = jj^/vphp, and 
where , A^ and A-^ are non-universal numerical coef- 
ficients. 



Appendix B: Spin-spin correlation functions 

Starting from the original spin representation (i.e., be- 
fore the rotation lA), the transverse spin-correlator is de- 
fined as 



C + - [T)^{TrSt+^ (T)5,r(0))^, 



(Bl) 
(B2) 



where the average is taken with respect to the total 
Hamiltonian of the system Eq. ([T]). The correlation 
function is a physical quantity that does not depend 
on our particular choice of representation. Using the 
transformation lA, the Hamiltonian is transformed as 
in Eqs. [25] and I26| while the spin operators become 
U^S+U = S'+e^-^'^^.-^o) (see Eq. 



241. Thus, 



(r) = (T. (U^e:^^U) [U^ St^JU) 



(B3) 



(B4) 



H 

(B5) 
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Imposing 7 = \f2, and near the strong-couplin g K ondo 
fixed-point, the Hamiltonian 'H, maps onto Eqs. (31 1 and 



( 32 1 (the conduction electron term Hp is not changed 
since it doesn't depend on spin operators). We therefore 
eliminate the spin-degrees of freedom in the correlation 
function and obtain 



(r) = (7Ve*^>+"Me-*¥''(o)\ 



-S'[v] 



(B6) 
(B7) 



where the averages are now taken with respect to the 
effective action Eq. ( 33 1 . 



1. Spin-correlators in the disordered Kondo-singlet 

phase 

In the disordered "Kondo" phase ^ | 1 Jh}' 
Heisenberg coupling in 5^ Eq. ( 35 1 is a suitable ex- 
pansion parameter to compute the correlation function. 
Explicitly 



= -■50 [ip] „iVi + >. {r)p-iv>i (0) 



/ V [ip] e-'^oM-'SffM 



/PM e-^o {e„.o nr=i 


J^dTj ^ (^^(e^'^'+i^^^)-*^'^^^) +h.c.) 


1 e^V.+r. 








) +h.c.) 


} 



42") dTidT2...dTn {^r^e*'^'+"(^)e-"^'+"(^")e''^'+"-i(^") 



g-i<,3i+ 1 (ti ) (ti ) g-ii^i (0) 



(B8) 
(B9) 
, (BIO) 



where we have truncated the perturbative expansion at 
leading order. This is to be expected, since the Heisen- 
berg term only couples nearest-neighbors, and therefore 
spins at a distance of n sites only become correlated at 
order ( J^) in the perturbative expansion. The product 

/7"^g»¥'i+n(-r)g-i¥'i+>i(-rn) ^ ^ _ Q-iVi+li-ri) giVii-ri) g~iVi{0)\ 

can be calculated using Wick's theorem. For com- 
pactness in notation, we define the (local) two-point 
correlation function 



FiT) = (Tre 



piV»n(Tl)g-iV™(T2) 



1 



(Bll) 

(at T = 0) , (B12) 



1 



where a is a numerical factor a — 4e "^""^ ~ 1.261..., and 
its Fourier transform 



Jo 



Eo 



exp 



Eo 



(B13) 



Then, the expression for C+ (r) compactly writes 

C+- (,) = / (#)" I e-'- [F , n>0 

\f{t). n = 

(B14) 

One particularly interesting case is the local dynamical 
correlation Cq (t) cx [cf. Eq. B12 , which encodes 



the properties of a local Fermi liquid. Another one is the 
static, non-local correlation 



c+- (0) 



1 e-"/«<= 



(B15) 



2 n+ 1 ' 

where we have defined the correlation length = 

'-(ft*: 



2. Spin-correlations in the ordered dissipative 
phase 

In the ordered phase, the dynamics of the spin-chain is 
effectively given by the action of the Goldstone mode Eqs. 
(47 1. With this Gaussian action, and using the saddle- 
point equations in Eq. (39 1 to express e*"^'^"^^ = -0* (t). 



we can calculate the spin-correlation in Eq. (B6l as 



C+- (r) = (jT,e*^'+"(^)e-'^'(°)) = ^r^,^-^^^^.ir) ^^Mo)^ ^ 

(B16) 



= (-!)"( 
= (-l)"e~5(^-[''.+"M-''.(o)l') 



-i^, + „ir)^i-ddO)^ 



(BIT) 
(B18) 



where the factor (—1)" comes from the antiferromagnetic 
correlations induced by the coupling > 0,1^ and where 



[A,+„ (r) - (0)]^ = 



1 - e 

— 2' 



ikao7i—iuJm'!' 



(B19) 
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We study this correlation function in two different limits: have 
the local limit n = , and r = 0. In the first case we 
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